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Background 



“Wake Ingestion Propulsion Benefits,” L. H. Smith, AIAA 
Journal of Propulsion and Power, 1993. 

Boundary Layer Ingestion (BLI) 

Propulsion has the potential for 
significant reduction (5-10%) 
in Aircraft Fuel Burn 

Previous system studies: 

Daggett, et al., NASA-CR-2003-212670 <s% 

Kawai, et al., NASA-CR-2006-2 14534 
Plas, et al., AIAA 2007-450 
Nickol, NASA-TM-2008-21 51 1 2 ,0%H 

Nickol and McCuller, AIAA 2009-931 Bum Reduction 

Recent system studies: 

Tillman, et al., AIAA invited pres., 2011 
Hardin, et al., AIAA-201 2-3993 


5% 



□ 



Boeing Boeing 

(Daggett, et al., 2003) (Kawai, et al., 2006) 


MIT NASA 

(Plas, et al., 2007) (Nickol, et al., 2009) 
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Technical Challenges 


• The potential benefits of Boundary Layer Ingestion (BLI) 
Propulsion can be diminished if key parameters do not 
meet their targets 

- Inlet total pressure loss 

- Fan efficiency reduction 

- Fan stall margin reduction 


- Fan aeromechanics requirements 
(dynamic stresses and flutter stability) 

v ) 
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Fan CFD Analysis - TURBO Code 


Implicit, finite-volume solver 

Reynolds-Averaged Navier Stokes equations 

Structured multi-block code 

Multi blade-row code 

k-epsilon turbulence model 

Inlet distortion boundary condition 

Throttle exit boundary condition 



Dynamic grid deformation for blade vibration 

Prescribed harmonic blade vibrations with energy method 
to evaluate flutter stability 
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Fan Computational Domain 


• Analysis of an Aero Design Iteration (not the Final Design) 




Fan rotor 
blade 


Exit 

boundary 


hub 


Inlet 


boundary 


Core-bypass 

splitter 


Exit 

boundary 


H, O, and C blocks 
of mesh generated by UTRC 
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Total pressure ratio 


Fan Performance - Clean Inflow 


• TURBO code (RANS solver) used with radial inlet profile 
of total pressure, total temperature, and flow angles 

• Speedline traversed by setting exit throttle condition and 
converging flow solutions 




Inlet corrected mass flow rate, lbm/s Inlet corrected mass flow rate, lbm/s 
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Inlet Flowfield Provides Distortion Pattern 


• Inlet flow computations were performed at UTRC for an 
inlet design iteration (not final design) and the flowfield 
results were provided to NASA 
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Fan Computation with Inlet Distortion 




Inlet distortion is prescribed as boundary condition at inlet 
boundary of the fan computational domain ( 18 -blade fan 
rotor and splitter) 
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Total pressure ratio 


Periodicity of Flowfield Around the Rotor 




Total pressure ratio for various blade passages 

1.55 
1.50 
1.45 
1.40 
1.35 
1.30 

29400 30400 31400 32400 33400 34400 


passage 1 
passage 5 
passage 10 
passage 15 


1 rotor 


revolution 



Time step counter 

\ 


Variation of total pressure ratio in different blade 
passages shows flowfield is converged to periodicity 




Total pressure ratio 


Periodicity of Flowfield Around the Rotor 


• Total pressure ratio for various blade passages 



/ 1 1 

Inlet Distortion causes variations in mass flow 

rate and pressure ratio around the fan rotor 


Structural Dynamics Model & Results 


Blade structural model created based on aero 
design iteration (structural design is in progress) 

• 8-node brick elements 

• 9,782 elements, 15,096 nodes 

• 222 nodes at the root constrained 



mode 1 mode 2 mode 3 mode 4 

63.5 Hz 156.6 Hz 224.8 Hz 346.6 Hz 
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Aeroelastic Formulation 



• Blade structural dynamics modal equations with 
aerodynamic load 

+[*]{<?} = { ad} 

aerodynamic load vector - 

Modal Force 



\ 

Modal Force computation 
requires unsteady pressure 
and modal displacements 





Forced Response 
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Modal Force 


Time history over one rotor revolution 




1,800 time steps 

per revolution 
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Modal force magnitude, lb 


Modal Force 


Fourier components 



Harmonic number 
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Campbell Diagram 
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Modal force magnitude, lb 


Modal Force 


Fourier components 
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Generalized displacement 


Forced Response - 

Vibrat i on A mplitude and Dynamic Stresses 



• Dynamic stresses are required to determine fatigue 
characteristics (Goodman diagram) 

{tfnr} = \[K n ] - 0i r 2 \M n ]1 { AD nr } forn th mode, r th harmonic 

dynamic stress o r = hs n q nr where s n is the modal stress 



harmonic or 
engine order 

vibration amplitude 
(inch) at tip t.e. 

dynamic stress 
amplitude (psi) 

1 

5.5 x 10- 2 

273 

2 

3.0x10- 2 

290 

3 

1.9x10-2 

666 

4 

3.1 x 10- 3 

308 

5 

2.6 xIO- 3 

169 

6 

2.7 xIO- 4 

33 

7 

7.0 x 10- 4 

427 

8 

6.0 xIO- 5 

19 
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Flow Chart for Flutter Stability Computation 



• Aerodynamic damping computation using TURBO-AE 




Mode Shape 



Calculate Aerodynamic 

Damping 
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Aerodynamic damping, % 


Flutter Stability with Clean Inflow 



Design operating speed, mode 1, 0 nodal diameter pattern 
(all blades in-phase), 18 blade passages (full rotor) 
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Aerodynamic damping, % 


Flutter Stability with Clean Inflow 


• Design operating speed, 18 blade passages (full rotor) 

• Phase angle of vibration = 360 * Nodal Diameter / 1 8 


mode 1 



Nodal diameter 
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Flutter Stability with Clean Inflow 







Design operating speed, 18 blade passages (full rotor) 
Phase angle of vibration = 360 * Nodal Diameter / 18 
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Aerodynamic damping, % 


Flutter Stability with Clean Inflow 


• Design operating speed, 18 blade passages (full rotor) 

• Phase angle of vibration = 360 * Nodal Diameter / 1 8 



Nodal diameter 
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Flutter Stability with Distorted Inflow 



Various Approaches 

• Circumferentially average the distorted inflow to obtain an 
equivalent radial profile; use work-per-cycle analysis 

• Select a portion of the inlet distortion to represent a “worst- 
case” inflow condition that is used at all circumferential 
locations; use work-per-cycle analysis 

• Prescribe blade vibrations and distorted inflow; use work- 
per-cycle analysis; average the results over all blades, and 
over multiple blade vibration cycles 

• Use tightly-coupled aeroelastic analysis with distorted 
inflow; blade vibrations are determined as part of the 
computations; post-process time history to estimate 
average damping over all blades and multiple vibration 
cycles 
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Flutter Stability with Distorted Inflow 



Current Preferred Approach 

• Prescribe blade vibrations and distorted inflow 

• Use work-per-cycle analysis 

• Average the results over all blades, and over multiple blade 
vibration cycles 


Work 


-<& f 

%/ rvrlp %/ < 


,d x 

-p.d A • ( )dt 

cycle surface ' — c)t 




Unsteady pressure includes effect of 
1 ) inlet distortion 


2) blade vibration c=> 


isolate this component to 
assess flutter stability 
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Aerodynamic damping, % 


Flutter Stability with Distorted Inflow 



• Design operating speed, mode 1 , 0 nodal diameter pattern 
(all blades in-phase), 18 blade passages (full rotor) 
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Vibration cycle number 


Blade number 


Aeromechanics Analysis of a Boundary Layer Ingesting Fan 


27 


Aerodynamic damping, % 


Flutter Stability with Distorted Inflow 



• Design operating speed, mode 1, 1 nodal diameter pattern 
(all blades in-phase), 18 blade passages (full rotor) 



Vibration cycle number 


Blade number 
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Summary 


• Created structural model based on aero design iteration 
and computed structural dynamics characteristics 

• Performed aeromechanical analysis of design iteration 

• Performed fan flutter analysis with clean inflow at design 
speed - no flutter encountered at conditions analyzed; 
additional work needed at part-speed conditions 

• Performed distorted inflow analysis for forced response 
vibrations to determine dynamic stress at design speed - 
additional work needed at on-resonance conditions near 
design speed 

• Performed initial analysis with blade vibrations and 
distorted inflow to estimate flutter stability - additional 
flutter analyses needed for other vibration modes and 
operating conditions 
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Future Work 


• Perform aeromechanical analysis on final inlet-fan design 
to ensure safe wind-tunnel test 

• Develop tightly-coupled aeroelastic analysis capability in 
TURBO for more detailed analysis of blade vibrations with 
distorted inflow 

• Perform aeromechanical analysis on updated fan stage 
design including non-axi-symmetric exit guide vanes 

• Develop inlet-fan coupled aeroelastic analysis capability 
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Questions? 
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